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Abstract 

We derive the lepton number violating dimension-five and dimension-seven operators, relevant 
• for neutrino mass generation, in the Minimal Supersymmetry Standard Model without i?-parity 

o . 

CN . (the -^MSSM) by using the effective Lagrangian method. We keep all the possible CP violating 

I phases, and we establish the general relationship between the high-scale parameters and the low- 

' energy observals associated with the Standard Model particles. We study in a specific model the 

dependence of neutrino masses on the CP phases. 
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I. INTRODUCTION 



Neutrino physics measurements show that neutrinos change their flavour in oscillationary 
manner Q|, implying that neutrinos have a mass. In the Standard Model (SM) of particle 
interactions neutrinos are strictly massless due to lepton number conservation and non- 
existence of right-handed neutrinos. The same is true for the most popular extension of the 
SM, the Minimal Supersymmetric Standard Model (MSSM) [2], where the lepton number 
conservation follows from the assumed R-parity symmetry. The R-parity of a field is defined 
as R = (^— l)3B+i+2s _ (^_]^~)3(B-L)+2s^ where B, L and s are the baryon number, lepton 
number and spin, respectively jsl. The conservation of R-parity is imposed in order to 
prohibit the proton from fast decay, which may happen if both B and L are broken. If one 
abandons the R-symmetry by allowing lepton number to be violated, but leaving baryon 
number unbroken, the MSSM can accommodate |AL| = 2 Majorana masses for neutrinos 



The superpotential of the MSSM with broken R-parity and conserved baryon number 
(acronymed as^MSSM) is in the most general case of the following form: 

W = W,^+ W^^ (1) 

with 

+K.u.L:Q.jDl . (2) 



The part W^^ is the superpotential of the standard R-conserving MSSM, whereas W^^ consists 
of renormalizable bilinear and trilinear couplings that violate the lepton number conservation 
and thereby the R-parity. In Q H^, are Higgs doublet superfields and Q, L are left- 
handed quark and lepton superfields, all transforming as SU{2) doublets, and U, D, E are 
right-handed quark and lepton superfields transforming as SU{2) singlets. The flavour of 
the lepton and quark fields is denoted by indices J, J, K, and the summation over gauge 
indices is implicit. The R-parity violating three-lepton coupling obeys A^^^ = — A^^^. 



The breaking of supersymmetry is due to interaction and mass terms of component fields. 
In addition to the soft terms present in the MSSM, one should now include also the terms 
that break lepton number. The most general form of the soft-term Lagrangian in our case 
is then 

= /:f + £f , (3) 



where 



+m H^H H^H, 



and 



+B^Hi, + {m%HlL^ + h.c.. (5) 

Here A^^ (a = 1, 2, ■ ■ • , 8), A^^ (a = 1, 2, 3), A^ denote the gaugino fields corresponding 
to the gauge symmetry S'f/(3) x SU{2) x f/(l) in an obvious manner. 

In this paper we shall study neutrino masses in the framework of the ^MSSM. The 
induction of neutrino masses by R-parity violation couplings have been extensively analyzed 
in the literature 0]. Nevertheless, in most of these earlier studies the possible CP-phase 
effects on the neutrino masses have been ignored. In our analysis, we will keep all relevant 
CP violation phases of the Lagrangian as free parameters, and we shall show that they do 
affect the neutrino masses. 

We will assume that the new physics, i.e. the physics that goes beyond the SM, is 
associated with an energy scale that is considerably larger than the electroweak scale of the 
SM, so that the new physics effects can be reliably analyzed by using the effective theory 
method. In that method, one integrates heavy fields out and keeps only the SM particles in 



the resulting effective theory [6|. The effects of the heavy fields are parameterized in terms 
of high-dimensional operators of the light degrees of freedom. Due to the lepton number 
violating terms in the full Lagrangian, there will appear dimension-odd operators in the 
effective theory [7] . After the electroweak symmetry breaking, these operators will induce a 
small (Majorana) mass for neutrinos, in accordance with experimental results j^. 

The paper is organized as follows. In Section |n] we will present preliminary considerations 
concerning Higgs field and fermion field doublets and define the SM Higgs doublet in terms of 
the original doublets and their vacuum expectation values. In Section lTTTl we derive, following 
the analysis of , the Wilson coefficients at the matching scale for the operators relevant for 
neutrino masses by integrating out the supersymmetric particles and the heavy Higgs doublet 
fields. The connection between the high-energy parameters and low-energy observables are 
also discussed in this Section. The results of numerical analysis and conclusions are presented 
in the last Section. 

II. PRELIMINARIES 

In the MSSM with no conserved lepton number there is no a priori distinction between 
the down-type Higgs boson and slepton doublets, which are assigned with identical quantum 
numbers. One can therefore freely rotate the weak eigenstate basis [Ha, Lj) (/ = 1, 2, 3) by 
an SU (4) transformation. The lepton number violating couplings depend on the basis one 
chooses. Nevertheless, the couplings of the SM Higgs field that appear in the low-energy 
theory should not depend on an specific choice of the basis. 

The Higgs boson and slepton doublets can be presented as follows: 




V 



H 



J 



d 
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V2 
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(/ = 1, 2, 3) 



(6) 



L: 
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Here v^, Vq, Vj denote the vacuum expected values (VEVs) of the neutral components of 
the doublets H^, and L^, respectively. One is obviously free to rotate these doublets in 
such a way that only one doublet achieves a non-vanishing vacuum expectation value. Let 
us denote 



(7) 



and replace the Higgs boson and slepton doublets of Eq. © with the set ($, 



$^^) defined as 



Hi 



(8) 



where denotes the transformation matrix 
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The doublets $, ^j:^^, (a = 1, 2, 3, 4) can be presented in the form 

( ] 
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(10) 



The doublet $ corresponds to the Higgs doublet of the SM, having a vacuum expectation 
value V — y/vf+vj = 246 GeV. It remains massless prior to the electroweak symmetry 
breaking and does not mix with the other doublets. The other doublets obtain a large mass 
originating in the soft supersymmetry breaking terms of the Lagrangian. Indeed, the masses 
of the doublets are given by 

/ $ \ 



1/ 



where 
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We have used here the following notations: 
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(12) 



(13) 



Here /, J = 1, 2, 3 and p,a — 0, 1, 2, 3, and we have denoted (m?)oo = mjj and 
(m^)jQ = (m^)o7. In the effective theory the heavy eigenstates of this mass Lagrangian are 



integrated out and only the fields in the massless doublet $ are considered as a physical 
degrees of freedom. 

The bilinear R-parity violating terms of the Lagrangian Q lead to mixing of the down- 
type higgsino and the left-handed leptons, and some of these fields obtain a large mass. The 
mass eigenstates are obtained by the following transformation: 

K I I (j^i^ 2, 3), (14) 




where the transformation matrix Z- is given by 



h 



- JfiLg-«'Pi 

1*^0^2 1 



l^ll 
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3I ^-i{Lp^-Lp^) 





Here the following notation is used: 



(15) 



\ P=0 



/^i^2 = 



I ^ l"=l I ) 



2 

3 1 ' 



(16) 



The phases (/ = 1, 2, 3), appearing in the transformation matrix, are defined as 
(yj^ = argie^ — 9^, where 9^ = arg{e^). The mass terms induced by the bilinear terms are 



(17) 



The phase 9 can be absorbed into the fields by the following redefinitions: 



(18) 



One can then deduce that Eq. p7p is a mass term of the Dirac field = yip„^ ' V'^ j • This 
heavy higgsino field, with a mass ^h, is integrated out in the effective theory. 
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FIG. 1: The tree level diagrams that lead to the dimension-five operators. 



III. THE EFFECTIVE LAGRANGIAN WITH DIMENSION-ODD OPERATORS 



In the effective theory that results when the heavy fields are integrated out neutrino 
masses are in leading order generated by the following dimension-five operators consisting 
of the SM fields Q: 

+liC,)UW^'^*){^^^'<) + h-c, (19) 

where (CJ^^, (Cj)^^ are Wilson coefficients. The tree-level Feynman diagrams inducing the 
above effective Lagrangian are given in Fig. ^ The Wilson coefficients corresponding to 
these diagrams are 

(c.),, = t^^o,K.^-''' > (20) 

where tan/3 = v^/v^, = cos/3, s„ = sin/3, and the parameters Aj j are definited as 



■e 



H 



^if^H 1=1 ^d/^ 

A^^^ = Jl^^Hl^ _ y ^ (21) 

with /, J = 1, 2, 3. When the EW symmetry breaks down, one ends up with the following 
mass matrix for neutrinos: 



X 



( JS? A A A A \ 

1 0,1 0,2 ■'^0,l-'^0,3 



■^0,1-^0,2 "^0,2 ■^0,2-^^0,3 



VAo,,Ao,3 Ao,2Ao,3 AJ^ y 
This matrix is diagonalized by the unitary matrix 

/ 



where 



V 



A A 



Ao.3 ^0,l^;.2 ^0,2 

A AA, A 



\,2 ^0.3 

A AA, A ' 



(22) 



(23) 



^= v/|Ao,J2 + |A |2 + |A 12 



A = VlAo,iP + |Ao,2P, 



(24) 



givmg 



1 o o/(7?e-^^i g^e-''^' 



xdiag(0, 0, A^) . 



(25) 



Hence only one neutrino flavor acquires a non-vanishing mass at tree level. This mass is 
given more explicitly by 



1 gle~'^^ gle-'^^ 



2^1 



\m^ I 



+ 



\m. 



X 



2 2 



V, 



oleol + E^JeJe*"^ 



1=1 



(26) 



By setting the phases zero and taking account of the approximations used, we find our result 
coincident with that of Ref. jlo| . 

An obvious conclusion from above is that the neutrino oscillation data, which indicates the 
existence of two distinctly different mass-difference scales, is not explained in the MSSM at 
tree level. Therefore, a study of higher order contributions to the neutrino mass Lagrangian 
is necessary. 



FIG. 2: The tree level diagrams that lead to the dimension-seven operators. 



Let us first consider the corrections of dimension-seven operators to the neutrino mass 
matrix. The tree-level Wilson coefficients of the dimension-seven operatos and the Wilson 
coefficients of the dimension-five operators at one-loop level belong to the same order of 
perturbative expansion. The dimension-seven operators that give corrections to the neutrino 
mass matrix are the following: 

+^(cj.,(<fV"<i.)((^i.*)(i.tav; 



2 

+h.c. , (27) 
where the Wilson coefficients have, refering to Fig. |21 the expressions 



(^5),/,/ — (Q),/,/ ' 



— I™ 12,, ^fl'^'p^ 



^KoKjKj ■ (28) 
After the electroweak symmetry breaking these operators yield corrections to the neutrino 
mass matrix, which are accounted for by the following replacement in Eq. (j22I): 

V I \m. I / 
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1 s r A 

^ 2 p 0,0 



+ (29) 

The same replacement should be done in Eq. (|26p. We realize that when the tree- level 
contributions of the dimension-seven operators taken account, still only one neutrino flavour 
acquires a nonzero mass. 

In Fig. 121 and Fig. 0] we present the one-loop diagrams that contribute to the Wilson 
coeflicients of dimension-flve operator after the matching procedure. In the full theory, 
the self-energy corrections to external leg, appearing in the diagrams in Fig. 01 should 
in princible include also the light-light contributions, such as e^e^, $$. However, such 
diagrams would not give any contribution to the dimension-flve operators after the matching 
of the full and effective theory jo | . A similar conclusion is true also for the vertex correction 
diagrams that involve only the light degrees of freedom. 

The lengthy expressions of the ensuing Wilson coefficients are presented in appendix 1X1 
After the EW symmetry breaking, the neutrino mass matrix with the one-loop corrections 
included can be written as 

M, + 5M,, (30) 

where denotes the contribution from the tree level dimension-flve and dimension- seven 
operators and bM^ is the correction from the dimension-flve operators at the one-loop level, 
given by 

+5Cf)) , (31) 



where the quantities bC^f^ are deflned in appendix 1X1 

The neutrino mass matrix is diagonalized, as usual, by a unitary transformation 

[/t (m^ +SM^)U = diag(m^^ , m,^ , mj , (32) 

where the 3x3 unitary matrix U generally contains 3 rotation angles and three physical 
phases (CP-phases), and denote the real positive mass eigenvalues of the light Majorana 
neutrinos. Among the three CP phases, two are so-called Majorana phases associated with 
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FIG. 3: The one-loop diagrams that lead to the dimension-five operators after the matching. The 
triangle, and self-energy diagrams denoted by blobs are presented in Figs. 5 to 8. 



lepton number violating terms in the mass Lagrangian and the third one is a so-called Dirac 
phase. One can factorize the Majorana phases by defining the transformation matrix U as 



L m ' 



(33) 



where is the lepton mixing matrix, a counter part of the Cabibbo-Kobayashi-Maskawa 
mixing matrix of quarks, containing the Dirac phase, and is a diagonal matrix with two 
Majorana phases. One can use the following parametrization: 

/ p^'^i ^ 
= e^'^a , (34) 
1 
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FIG. 4: The box-diagrams that lead to the dimension- five operators after the matching. 



A5 



j<5 \ 



(35) 



V 

with c = cos 6*.., s.. = sin 6'.. and 9.. denote the mixine angles. 

The rotation angles 6ij and the phases 6 and Ui appearing in the lepton mixing matrix 
Ul, as well as the values of neutrino masses, can be expressed in terms of the parameters 
of the the model. In the following we will study the constraints set by the atmospheric and 
solar neutrino oscillation data on the model parameters. 



IV. CONSTRAINTS FROM OSCILLATION DATA 



The neutrino oscillation data have outlined the mass and mixing pattern of neutrinos. 



According to a global fit to t 
following set of parameters ^| 



le data, the oscillations are best explained in terms of the 



Aml^ = 2.0 X 10"^ eV^ 



Am^ = 7.2 X 10^^ eV^ 

21 
2 , 



sin 6*23 = 0.5, sin 9-^^ = 0.3, 
sin2^,3 < 0.06 . (36) 



Let us now discuss the implications of these empirical results for the parameters of the 
-^MSSM. Beside the parameters appearing in the MSSM, the -^MSSM includes generally 
many new ones (in general complex) related to lepton number violation: three bilinear (e^) 
and thirty six trilinear (nine A^-^^, twenty seven A'^^^) couplings in the superpotential, to- 
gether with the corresponding soft breaking (three B^) and A-terms (nine A^j^, twenty seven 
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FIG. 5: The one-loop diagrams that lead to the dimension-five operators after the matching. The 
triangle-, and self energy-diagrams are presented in Figs. 5 to 7. 

A'^^^) and three additional-^ soft masses (m^),,^. There are only six physical parameters 
among the nine bilinear parameters (e^, B^, (mpg^) due to the ambiguity in the choice of 
the {H^, L) basis. It is obvious, giving the large number of parameters, the oscillation data 
alone would leave in the general case much free space for the choice of parameters. In the 
following we will not consider the most general situation but restrict ourselves to the case 
where only the bihnear.^ couplings are present. We are particularly interested in the effects 
of the complex phases associated with these couplings. 

Neutrino masses in the framework of.^MSSM are extensively discussed in the literature, 
but the studies are usually restricted to the CP-conserving case leaving the possible effects 
of CP phases with less attention. In the MSSM, the soft broken terms provide new sources 
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FIG. 6: The gaugino self-energy diagrams that lead to the dimension-five operators after the 
matching. 

for the CP violation, in addition to the CKM mechanism of the Standard Model. At present, 
the strictest constraints on the CP phases originate from the experimental bounds on the 
electric dipole moments (EDM) of the electron and the neutron. Nevertheless, if one invokes 

n 

a cancellation mechanism among different supersymmetric contributions [Ij] or choose the 
sfermions of two generations heavy enough the loop- induced EDM's yield a bound only 
for the argument of the parameter n, implying /j, < 7r/(5tan/5), leaving the other explicit 
CP- violation phases unconstrained. 

In principle, all the parameters associated with i?-parity violation are complex in the 
model we are considering. For simplicity, we will now assume that all parameters are real 
except the gaugino masses. Furthermore, in order to apply the effective theory method 
safely, we assume the new physics scale is fi^p = 10 TeV, in other words, much higher than 
the electroweak scale. We choose the basis t;^ = (/ = 1, 2, 3), and we set 

(m-),, = diag(l, V2, 73) x 10 TeV, 
(m-)„, = (100, 100, 100) GeV, 

i^E^ij = i^u^iJ = (^D^iJ = diag(5, 5, 1) TeV, 

= diag(5, 5, 1) TeV, 
Al = M=Al = 0,{I, J=l, 2, 3). 
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FIG. 7: The lepton higgsino self-energy diagrams that lead to the dimension-five operators after 
the matching. 





FIG. 8: The Higgs self-energy diagrams that lead to the dimension-five operators after the match- 
ing. 
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For the bilinear couplings we use the values 



= 5 TeV, = 0.01 GeV, 
= 0.3 GeV, = A GeV, 
= 100 TeV^ B^ = B.^=B^ = 100 GeV^ . 

(37) 

In order to obtain small neutrino masses, < 0.1 eV, we assume the gaugino masses very 
large, |mj = Im^l = 100 TeV. We consider the numerical values of the parameters given 
above to be conceivable and representative for the model. 

In order to demonstrate the effects of CP-phases on the masses of neutrinos we plot in 
Figini in the case 9^ = arg{m^) = 0, the neutrino mass-squared differences Am^^, Am^^, 
and the neutrino masses fn^^, m^^ versus the CP-phase 6-^ = arg{m^). We have considered 
the situation for both large (tan/3 = 50, solid line) and small (tan/3 = 5, dashed line) 
values of the ratio tan/3 = Vu/vd- In FigEl we present a similar plot the same for 9^ = 
arg{m^) = and a varying 6^ = arg{m^). It is seen from these plots that the solar and 
atmospheric neutrino data, which allow the 3o"-range shown in gray in the figures, sets 
non-trivial connstraints on the complex phases, in particular in the case of a small tan (3. 

Let us finally note that we have found, by scanning the parameter space, that when all 
the trilinear ^ couplings together with corresponding soft terms are included, it is possible 
to realize in this model a situation, where a large neutrino mass hierarchy is associated with 
a large mixing. In general, this kind of situation is not so easy to achieve in a definite 
theoretical framework. 

V. SUMMARY 

There are many ways to generate neutrino masses in the Minimal Supersymmetric 
Standard Model without R-parity. In this paper we have derived, using the effective field 
theory approach, the high- dimension operators relevant for neutrino masses, which the R- 
parity breaking couplings induce. Among these high-dimension operators, the dimension- 
odd operators originating from the lepton number violation couplings will lead to nonzero 
Majorana masses after the EW symmetry breaking. By properly defining the heavy fields, 
we have derived the Wilson coefficients of dimension-five operators at one loop level and the 
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FIG. 9: The neutrino mass-squared differences Am?,, Am?,, and the neutrino masses m„ , m„ 
versus the CP phase = arg{m^) for O,^ = arg{m^) = 0. The sohd hnes correspond tan/3 = 50 
and the dash lines to tan /? = 5. The gray bands represent the region allowed by the global fit of 
the solar and atmospheric neutrino data at Scj level. For the values of the other parameters, see 
the text. 

those of dimension-seven operators at tree level. We have demonstrated the effect of CP 
phases on neutrino masses by looking at a specific case, where the only complex i?-parity 
violating couplings are the soft bilinear gaugino mass terms. 
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FIG. 10: The neutrino mass-squared differences Amg^, ^'^ai' ^^'^ neutrino masses "m^^, m^^ 
versus tlie CP pliase 6^ = argim^) for 9^ = arg{m-^) = 0. The sohd lines correspond to tan j3 = 50, 
the dash hues tan /J = 5. 



APPENDIX A: THE ONE-LOOP CORRECTIONS ON THE WILSON COEFFI- 
CIENTS OF THE DIMENSION-FIVE OPERATORS 

In this appendix, we present the one-loop corrections on the Wilson coefficients of the 
dimension-five operators that are present in Eq. (jl9j) . The triangle diagram corrections are 



2 m, 
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with X. = m^/u^,„. 

» i I r~ NP 

The corrections from the self-energy diagrams 
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where x. — / /i'^^ , and the nonzero Yukawa couphngs are 
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Sii )/,2 
Sii )l,3 



h'^ leJe^'^i - A' leJ 



Hi 



ff2 



with 



V, 



pij 



The couphngs between the Higgs and right-handed slepton are given by 

K-l 



(A3) 



(A4) 



with 



SjJO,a 



3 3 ^1 _ vk -'^K y^vi 



K=l 



ii"=i 



vk 



J=l 



V 
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+ 



0.2 a 



V. V. V. V, V, V. 



V 



d2 



The loop-integral functions are 

1 — Inx 



A{x) = 



Ao{x) = X 
B{x,y) = 
Bo{x,y) = 
P{x,y) 



(47r)2 ' 
1 — Inx 



(47r)2 ' 



1 xlnx ulny 
1 + + ^ — - 



(47r)2 
1 

(4^ 



y — X X — y 

„2 



a; Inx y Iny 



y — a; a; — 7/ 



1 p Inx ^ Iny 



(47r)2 ^x — y y — X 
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C{x,y,z) 



1 r a; In a: 



(47r)2 '-{y — x){z — x) 



{x-y){z-y) {x-z){y-z) 
1 r x'^lnx 



\ny z\nz 



Co{x,y,z) 




{x-y){z-y) {x-z){y-z) 



- 1 
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